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ABSTRACT 

We propose that generic magnetic equilibrium of an ideally conducting fluid contains a volume- 
filling set of singular current layers. 

Singular current layers should exist inside neutron stars. Residual dissipation in the singular 
current layers might be the main mechanism for the magnetic field decay. The slow decay of the 
field might be the clock responsible for triggering the magnetar flares. 



1. The magnetic equilibrium problem 

An equilibrium (no convection) star might still 
be a good conductor and contain magnetic fields 
- with the magnetic stress balanced by the fluid 
pressure and gravity. We want to describe the 
generic features of such equilibrium magnetic 
fields. 

This problem was clearly formulated many 
years ago by Chandrasekhar and Fermi, but it 
remains unsolved to this day. Recent progress 
came from numerical work (Braithwaite & Nord- 
lund 2005). 

The recent numerical work does not point out 
an important generic feature of the equilibrium - 
singular current layers. Moreover, we propose that 
generic equilibrium carries infinitely many singu- 
lar current layers. 

This proposition - Generic equilibrium carries 
a volume-filling set of singular current layers - is 
our main thesis. This is indeed just a proposition. 
It is based on very shaky arguments. A numerical 
proof or refutation is needed. 

The author tried to solve the problem numeri- 
cally, but succeeded only in 2D, where the singular 
current lines are isolated rather than area-filling. 
Still, it appears that modern computers, if han- 
dled by professionals, are powerful enough to also 
solve the 3D problem (by the relaxation technique, 
§2). 

Since it is believed that close to equilibrium 
magnetic fields do exist in neutron stars and in 
other non-convective stars, a careful numerical 



study of the singularities should be worthwhile. 
If singularities are indeed generic, they might be 
responsible for some of the observed phenomena - 
say, magnetic field decay and triggering the mag- 
netar flares. 

The singular current layer is a local feature, it 
can be considered without gravity. To simplify 
the analysis, we will also assume an incompressible 
constant density fluid. Then the magnetic fields B 
will be in equilibrium if the magnetic j x B force 
is balanced by the pressure gardient: 

B X V X B = -Vj9. (1) 

We merely want to discuss solutions of eq.([l]). 

2. The relaxation method 

A well-known procedure for finding stable mag- 
netic equilibria is the relaxation method. Con- 
sider a viscous but ideally conducting incompress- 
ible fluid. Set an arbitrary initial magnetic field 
and velocity, and let the system evolve. The total 
energy will decrease due to internal friction, and 
asymptotically the fluid must stop. The fluid can 
stop only in a magnetic equilibrium state. Since 
the fluid is an ideal conductor, the magnetic field 
topology is not changed by the fluid motion. Con- 
served topology plus incompressibility guarantee 
that, generically, the equilibrium state reached by 
the fluid will be non-trivial, B 7^ 0. 

The relaxation method proves the existence of 
magnetic equilibria with arbitrary topology. The 
method is constructive - it allows one to actually 
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find this equilibrium. 

The corresponding mathematical problem, in 
the limit of large viscosity, and after an obvious 
generalization reads: 

B = V X (v X B), 
V = -LV X V X (B X V X B), ^ ^ 

where L is an arbitrary positive operator which 
commutes with V. The first equation says that 
the magnetic field B is frozen into the fluid moving 
with velocity v. The second equation says that the 
fluid incompressibly yields to the magnetic force; 
L = corresponds to internal friction, L = 

—A~-^ corresponds to external friction (A = V^). 

Eq.([2]) describes a dissipative evolution. Mag- 
netic energy 

W = ^J d^rB^ (3) 

decreases: 

= - / (V X (B X V X B)) 

•L(V X (B X V X B)) < 0. ^ ^ 

The evolution stops only when 

V X (B X V X B) = 0, (5) 

which is equivalent to eq.([l]). 

We must conclude that for any solenoidal field, 
there exists an incompressible deformation of this 
field which is an equilibrium. 

3. The Arnold argument 

The following argument, due to Arnold (1986), 
explains why the magnetic equilibrium problem re- 
mains unsolved, despite the clear and simple logic 
of the relaxation method. 

From eq.([l]), B • Vp = 0. This means that mag- 
netic field lines lie on the surfaces p=const. But 
the magnetic field we started with in the relaxation 
method generically does not possess magnetic sur- 
faces - generic field lines fill up volumes, not sur- 
faces. Clearly, the isotopological relaxation ([2]) 
cannot comb the intertwined volume-filling lines 
onto the nested tori p=const. 

And yet, the relaxation argument of §2 is still 
valid. For any solenoidal field there must ex- 
ist at least one incompressible deformation of the 



field which is an equilibrium, because asymptoti- 
cally, the motion described by eq.([2]) stops. It also 
appears likely that for suitable operators L, the 
relaxing field remains non-singular at any finite 
time, meaning that the relaxation method indeed 
stops at a well-defined asymptotic magnetic field 
(the final field may be singular and different for 
different operators L). 

4. The 2D relaxation and equilibria 

After the above discussion, one really wants to 
offer this problem to a computer and see how the 
computer will resolve this contradiction between 
topology and equilibrium. The author was able to 
obtain trustworthy numerical results only in 2D. 
In fact, the 2D case is reach enough to show an 
important generic feature of magnetic equilibrium 
- singular current layers. We discuss the 2D equi- 
libria in this section. 

First of all, it appears that in 2D there are 
no topological obstacles to equilibrium. Indeed, 
Arnold's argument says that magnetic field lines 
lie on the lines of constant pressure. But the 2D 
field lines are always the isolines of the magnetic 
stream function: V • B gives B = {—dyip^ dxi^)- 

Therefore one expects that a smooth equilib- 
rium state will correspond to any initial state ip. In 
reality, it has been known for some time that ini- 
tial states with non-trivial topology relax to equi- 
librium states with singular current layers 

To show how the singular current layers form, 
consider relaxation with internal friction. L = 



The author does not know who discovered that singular 
current layers are generic; would be grateful for a reference. 




Fig. 1. — Left: Initial isolines of tp. Right: Final 
isolines of tjj. 
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Fig. 2. — Left: Initial isolines of tjj. Right: Final 
isolines of t/^, also shown (thick colored lines) are 
the isolines of At/j. 




Fig. 3. — Graphical "proof that singular current 
layers form from the X-points 




Fig. 4. — Left: Initial isolines of ip. Right: Final 
isolines of tp. Also shown (thick colored line) is 
the theoretical separatrix - two 251° circular arcs 
and a line segment. 



A ^ in eq.([2]) gives the "incompressible diffusion" 
equation in 2D: 

{/, g} = da:fdyg - dyfd^g. 
The magnetic energy decreases 

W=^l d^r{Wi^)\ W = -J dV(Ax)^ (7) 

The evolution stops at equilibrium 

{^,A^} = 0, (8) 
or, in another form, 

All; = F{tlj). (9) 

Computer simulations of eq.(|6]) (using fast 
Fourier transform) show that all initial saddles 
of i/j develop into singular current layers, see fig- 
ures (1) and (2). The singularities are of the type 
i/j (X giving the field B ex sign(x), and the 
current j (x S{x) - singular current layers. 

The reason for existence of singular current lay- 
ers in 2D appears to be clear. The incompressible 
relaxation process (|6]) tries to make the field lines 
as short as possible, while keeping the area en- 
closed by each field line constant. The field lines 
try to deform into circles. This is impossible for 
non-trivial initial topology. Then, as illustrated 
by fig.®, one does expect that the saddles (X- 
points) of the initial field split into two Y-points 
with a singular current layer between them. H 

One can actually transform the cartoon fig.® 
into a rigorous calculation. Take the limit in which 
only an infinitesimal layer of the field surrounds 
the separatrix. Assume equal initial energies in 
the outer and inner parts of the separatrix. Then, 
the magnetic energy will scale as W ex L| + L^, 
where Li is the inner and Lq is the outer lengths 
of the separatrix. Minimization of the energy at 
fixed enclosed areas predicts the final shape of the 
separatrix - two ~ 251° circular arcs and a line 
segment separating the arcs, for a symmetric sep- 
aratrix. The direct numerical simulation confirms 
this result, fig.®. 



^ Does the outer part of the separatrix become a singular 
current layer? In the past, the author claimed that it does 
but the argument given was incorrect. Our numerical re- 
sults are inconclusive. 
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5. Perturbation theory in 2D and in 3D 

We have also performed 3D numerical simula- 
tions of the incompressible relaxation ([2|) with the 
inner friction, L = A~^. We see clear evidence for 
the singular current layers formation!! but we were 
unable to see how many singular layers form, how 
are they distributed in the volume, and how ex- 
actly does the system manage to preserve the field 
line topology, while putting the field lines onto the 
constant pressure surfaces. One probably needs a 
better than our 64^ resolution, and/or a better 
handling of the numerical issues. 

Having failed numerically, we resort to linear 
perturbation theory. We start with a translation 
invariant equilibrium field, put a small perturba- 
tion on it, and ask what does the perturbed field 
relax into. 

While this linear problem can indeed be solved 
analytically, it is unclear what does the solution 
have to do with the generic equilibrium, which, of 
course, is fully nonlinear. We will see, however, 
that in 2D the linear perturbation theory can be 
used as a sensible guide to the properties of the 
generic equilibrium. We will assume that the same 
is true in 3D, and this will allow us to propose that 
the generic magnetic equilibrium carries a volume- 
filling set of singular current layers. This kind of 
reasoning is unsatisfactory. But this all we can 
offer. 

5.1. Perturbation theory in 2D 

In 2D, the initial y-translation invariant back- 
ground field is ^po{x)^ corresponding to the mag- 
netic field B{x) = i^oix) in the y-direction. The 
perturbation is ex e*^^, and we assume that k 
- otherwise the perturbation simply re-difines the 
background field V^o(^)- The final relaxed state 
will be ipo{x) -\-il){x)e^^y . To linear order, the equi- 
librium condition ([8]) gives 

B^P" - B"^P - k'^Bijj = 0. (10) 

If B everywhere, one can write the perturba- 
tion 8iS ip = B(j). Then 

{B^^'y - k^B^ = 0. (11) 



^ Numerical 3D singular current layers have been known for 
many years. 



With the periodic or decreasing x-boundary con- 
ditions, this gives 

J dxB^{(^'^ ^k^(f^) =0, (12) 

and = 0. This result is trivial - just says that 
the field lines straighten up during the relaxation. 

Now consider the case when the background 
magnetic field B has zeros. In the short wave- 
length limit, for B{xo) = 0, eq. (p!Q|) reads 

{x - xo){tlj'' - k^i/j) =0. (13) 

The solution is 

i/j" - k'^t/j (X S{x - xo), . . 

This is a current singularity. 

So far everything was fully rigorous. What is 
shaky, is our extrapolation of this result to the 
non-linear regime. Based on the linear perturba- 
tion theory, we "conclude" that the generic 2D 
equilibrium will have isolated singular layers at 
magnetic field zeros. This is indeed correct in 2D, 
if by zeros we mean only the saddle-point zeros. 

5.2. Perturbation theory in 3D 

In 3D, we choose a background field of the 
form Bo = (0, /7(x), V(x)). The perturbation is 
= (b^.by.b^) (X e*^^+^^^ Then the linearized 
equilibrium condition ([T]) is 

i{mU + kV)b:, - {Uby + Vb,y = p^ (15) 

iV{kby — mbz) + U'bx = imp, (16) 

- iU{kby - mbz) + Vb^ = ikp. (17) 

Here ' = d/dx, and p is the pressure perturbation. 
The magnetic field perturbation is solenoidal 

b'^ + imby + ikbz = 0. (18) 

From (p!6|) and (p!8|) we calculate the product 
Vbz in terms of bx and p. From (pT|) and (p!8|) we 
calculate the product Uby in terms of bx and p. 
Plug these into ([TSj). The pressure perturbation p 
cancels, and we get 

Wb'^ - W"bx - (m2 + e)Wbx = 0, . . 
W{x) = mU{x) + kV{x) ^ ^ 
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This equation is isomorphic to eq.([TO]), and we can 
use the 2D results. If W{x) has no zeros, the per- 
turbation is zero. For the relaxation problem it 
means that the lines straighten up. But if W does 
have zeros, each zero will develop into a singular 
layer. 

The difference between 2D and 3D is that in 
3D the zeros of W{x) depend on wavenumbers. 
Suppose we work in a 27r-periodic cube. Then the 
wavenumbers m, k are arbitrary integers, and all 
points X where U {x) /V{x) is a rational number are 
zeros of W for the appropriate choice of m and k. 

"Based" on the above, we propose that generic 
magnetic equilibrium contains a volume-filling 
(meaning everywhere dense within some domains) 
set of singular current layers. 

This proposition leaves many questions unan- 
swered: 

• Is this proposition true at all? 
And assuming it is true: 

• How exactly does this infinite set of singular 
layers help resolve the topological obstacles 
to relaxation? 

• Where do the layers form? 

• Which current layers, a few strongest ones, 
or a lot of weaker ones, are responsible for 
the residual damping, once a small (probably 
anomalous in the real world) resistivity is 
added to the problem? 

6. Summary 

We propose that generic magnetic equilibrium 
contains a volume-filling set of singular current 
layers, because topological obstacles to relaxation 
are everywhere dense. This is just a proposi- 
tion. The proposition leaves many questions unan- 
swered. The proposition might turn out to be 
completely false. 

It would be interesting to see what a generic 
field relaxes into in 3D. 
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^ In the perturbative calculation, the singular layers form at 
rational U/V, that is on closed field lines. 
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